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Abstract 

We present measurements of the meson total cross section and differential cross section 
da/dpT- The measurements use a 98 ± 4 pb~^ sample of pp coUisions at ^/s — 1.8 TeV 
collected by the CDF detector. Charged B meson candidates are reconstructed through the 
decay B"^ — > J/t/jK^ with J/z/' — > A*"*"/^"- The total cross section, measured in the central 
rapidity region \y\ < 1.0 for Pt{B) > 6.0 GeV/c, is 3.6 ± 0.6(stat ® syst)/xb. The measured 
differential cross section is substantially larger than typical QCD predictions calculated to 
ncxt-to-lcading order. 

PACS numbers: 13.85.Ni, 12.38.Qk, 13.25.Hw, 14.40.Nd 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) can be used to compute the expected cross sections 
for the production of heavy quarks at hadron coUider energies. Calculations of the hard- 
scattering cross section have been carried out to next-to-leading order in perturbation the- 
ory [0. Experimental measurements must show that these predictions provide an adequate 
description of the cross section at 1.8 TeV before they can be confidently extrapolated to 
higher energies or more exotic phenomena. Unfortunately the QCD predictions are affected 
by large theoretical uncertainties such as the dependence on the choice of the factorization 
and renormalization scales, the parton density parameterization and the b quark mass 0. 

Experiments at CERN and at the Tevatron |Q have shown that the b quark production 
cross section is higher than the theoretical predictions obtained with the standard choice 
of parameters by about a factor of 2-3. Closer agreement between theory and the experi- 
mental measurements can be achieved by choosing rather extreme values of the theoretical 
parameters 0. It has also been suggested that the large discrepancy could be explained by 
pair production of light gluinos that decay into bottom quarks and bottom squarks 0. 

This paper describes a measurement of the meson total cross section and differential 
cross section da/dpr in hadronic collisions using fully reconstructed mesons decaying 
into the exclusive final state J/ip K^. The measurement uses a data sample of 98 ± 4 pb~^ 
collected by the Collider Detector at Fermilab (CDF) experiment from pp collisions with 
a center-of-mass energy of 1.8 TeV produced by the Fermilab Tevatron. The data were 
collected in the run period from 1992 to 1995 which is referred to as Run 1. Our previously 
pubhshed result based upon 19.3 ± 0.7 pb~^ of data (Run lA) found that the total cross 
section for priB) > 6.0 GeV/c and \y\ < 1.0 is as = 2.39 ± 0.54(stat © syst) /ib. 

The paper is organized as follows. In Section II we review previous measurements of the 
B cross section using exclusive B decays. In Section III we briefiy describe the components 
of the CDF detector relevant to the analysis presented in this paper. The data collection, 
event selection procedures and the reconstruction of 5^ J/ ipK^ are discussed in Section 
IV. The measurement of the differential and total cross sections is presented in Sections [V| 
and respectively. 
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II. PREVIOUS MEASUREMENT OF THE B PRODUCTION CROSS SECTION 

The Run lA measurement of the B meson differential cross section was determined 
from fully reconstructing the decays 5± J/tpK^ and 5° J/tljK*^{892) The 
measurement of the B transverse momentum spectrum showed that next-to-leading-order 
QCD adequately described the shape of this distribution for pr > 6.0 GeV/c. In the Run 
lA publication, CDF used a branching ratio BR{B+ J/i)K+) = (11.0 ± 1.7) x 10"'' and 
a product of branching fractions B = BR{B^ J /ipK^) x BR{J /ip = (6.55 ± 

1.01) X 10~5 [g. The current world average for BR{B+ J/ijK+) is (10.0 ± 1.0) x 10"^ 
which yields B = (5.88 ± 0.60) x 10"^. The change in the branching fractions scales the 
published result up by about 10% to asipx > 6.0 GeV/c, \y\ < 1.0) = 2.66 ± 0.61 (stat © 
syst)yub. 

This paper updates the measurement presented in 1995 by using the complete Run 1 
data sample of 98 ± 4 pb~^. For this measurement, we use only the decay mode B^ — >■ 
J/ipK^ where we require both muon candidates from the J/ip decay to be well measured 
by the silicon vertex detector (SVX). Such a restriction allows us to use fewer selection 
requirements since the decay mode B^ J/ipK^ has a lower combinatorial background 
than B^ — > J/ipK*^, and the SVX information enables us to substantially reduce the prompt 
background. Moreover, several of the efficiencies are measured using a large sample of J ftp 
^~^fi' candidates rather than relying on Monte Carlo calculations for detailed modeling of 
detector effects. 

III. THE CDF DETECTOR 

The CDF detector is described in detail in P). We summarize here the features of the 
detector subsystems that are important for this analysis. The CDF coordinate system has 
the z axis pointing along the proton beam momentum, and the angle is measured from 
the plane of the Tevatron storage ring. The transverse (r-0) plane is normal to the proton 
beam. 

The CDF experiment uses three separate detectors for tracking charged particles: the 
silicon vertex detector (SVX), the vertex detector (VTX), and the central tracking chamber 
(CTC). These devices are immersed in a magnetic field of 1.4 Tesla pointed along the —z 
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axis generated by a superconducting solenoid of length 4.8 m and radius 1.5 m. 

The innermost device is the SVX which provides spatial measurements in the r-cj) 
plane. The SVX consists of two cylindrical barrels that cover a region 51 cm long in z. Each 
barrel consists of four layers of silicon strip sensors with strips oriented parallel to the beam 
axis. The distribution of the pp collisions along the beamline is Gaussian in z with a a of 
about 30 cm. Therefore only about 60% of all J/ip ^ events have both muon tracks 

reconstructed in the SVX. 

The SVX is surrounded by the VTX, a set of time projection chambers which measure 
the z coordinate of the pp interaction (primary vertex) . Surrounding the SVX and the VTX 
is the CTC The CTC is a 3.2 m long cylindrical drift chamber with 84 layers of sense 
wires ranging in radius from 31 cm to 133 cm. The combined momentum resolution of the 
tracking chambers is Spr/pT = [(0.0009 + (0.0066)^]^/^ where pr is the component of the 
momentum transverse to the z axis and is measured in GeV/c. Charged track trajectories 
reconstructed in the CTC that are matched to strip clusters in the SVX have an impact 



parameter resolution of adipr) = (13 + 40/pt) /xm |]10| with px in units of GeV/c. The track 
impact parameter d is defined as the distance of closest approach of the track helix to the 
beam axis measured in the plane perpendicular to the beam. 

The central muon system consists of three components (CMU, CMP and CMX) and 
detects muons with p^ > 1.4 GeV/c in the pseudorapidity range |?7| < 1.0. The CMU 
system covers the region |?7| < 0.6 and consists of four layers of drift chambers outside the 
hadron calorimeter. Outside the CMU there is an additional absorber of 60 cm of steel 
followed by four layers of drift chambers (CMP). The CMX system extends the coverage to 
pseudorapidity 0.6 < |?7| < 1.0 but is not used in this analysis. 

CDF employs a three level trigger system. The first two levels are implemented in custom 
electronics. To select events in the third level, we employ a CPU farm using a version of the 
CDF event reconstruction program optimized for speed. 



9 



IV. DATA SAMPLE SELECTION 



A. Dimuon Trigger 

The data sample consists of events that pass the J/V' — > fi'^fi~ trigger. In the first level 
of this trigger, we require two muon track segments in the central muon chambers separated 
by at least 5° in azimuth. The trigger efficiency for each muon at Level 1 rises from 50% for 
Pt = 1.7 GeV/c to 95% for pr = 3.3 GeV/c. 

In the second level, we require muon segments found in Level 1 to be associated with 
tracks identified by the Central Fast Tracker (CFT) []rT|. The resolution of the CFT is 
Spx/Px ~ 0.03 (GeV/c)~^ In Run lA and for a subset of the Run IB data, we required one 
of the two muons to be matched to a CFT track with px greater than about 3 GeV/c while 
in the bulk of the Run IB sample, we required two muon segments to have an associated 
track with a threshold of about 2 GeV/c. In Run lA (IB), the extrapolation of the track 
was required to be typically within 10° (5°) of the muon segment. The efficiency of the track 
requirements was measured in a J/ip data sample using events in which the muon under 
study need not have satisfied the requirements for the event to be accepted. The efficiency 
for the nominal 2 (3) GeV/c threshold rose from 50% of the plateau efficiency at 1.95 (3.05) 
GeV/c to 95% of the plateau efficiency at 2.2 (3.4) GeV/c. That plateau efficiency changed 
over the course of the run because of aging of the CTC and subsequent modifications to the 
CFT algorithms. That dependence on time is accounted for in the calculation of the trigger 
efficiencies. 

The Level 3 software trigger required two muon candidates with an effective mass in the 
J/ip mass region after full reconstruction. Runs with known hardware problems for muons 
were removed yielding for this analysis a total Run 1 luminosity of 98 pb~^. 

B. J ftp Reconstruction 

Background events in the dimuon sample collected with these triggers are suppressed 
by applying additional muon selection cuts. Track quality requirements are used to reduce 
the backgrounds arising from poor track measurements. Tighter cuts are imposed on the 
correlation between the track in the muon chamber and the extrapolated CTC track. 

The transverse momentum of each muon from the J/tp for Run lA is required to be 
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greater than 1.8 GeV/c with one muon of the pair greater than 2.8 GeV/c. For Run IB, 
both muons are required to have a transverse momentum greater than 2.0 GeV/c. Events 
passing both the trigger and px requirements identical to those of Run lA are also accepted. 
The muons must have opposite charge and the separation in z between the two tracks must 
be less than 5.0 cm at the point of closest approach to the beamline. The z coordinate of 
the decay vertex is required to be within ±60 cm of the detector center. 

The invariant mass and uncertainty {(Tm) of the J/tp candidates are calculated after 
constraining the two muon tracks to come from a common point in space (vertex constraint) 
to improve the mass resolution. The width of the reconstructed J/ ip mass peak is 16 MeV/ c^. 
The signal region is defined to be those dimuon candidates with reconstructed mass within 
3.3am of the known J/tp mass |T^. We find (8.7 ± 0.2) x 10^ J/tp over background. In this 
analysis, the two muons from the J/ tp decay are required to be reconstructed in the silicon 
detector. 



C. Primary Vertex Selection 

Knowledge of the distance between the primary pp interaction vertex and the secondary 
decay vertex in the transverse plane is crucial to this analysis since the B meson proper 
lifetime is used to discriminate between B mesons and background events. We find the 
transverse position of the primary vertex using the average beamline calculated for each 



Tevatron pp store [T^. The longitudinal coordinate of the primary vertex (z) is measured 
using data from the VTX detector. The slopes and intercepts of the run-averaged beam 
position are combined with the event-by-event z locations of the vertices to determine the 
vertex position. The primary vertex uncertainties a^^, (Jy and are estimated to be 25, 25 
and 300 /xm, respectively. 



D. B Reconstruction 



To select charged B candidates we considered each charged particle track as a kaon 
candidate to be combined with a J/ip. A charged track in an event is combined with the 
two muons if the zq parameter of the track is within 5 cm of the z position of the J/ip 
candidate decay vertex. The exit radius of the kaon candidate, which corresponds to the 
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radius at which the track trajectory intersects the plane of the CTC endplate, is required 
to be greater than 110 cm to hmit the search to a region of high tracking efficiency. A 
cut on the kaon transverse momentum of Pt > 1-25 GeV/c is imposed to reduce the large 
combinatorial background. This cut is very effective since kaons from B meson decay have a 
considerably harder pt spectrum than particles from the underlying event and from events 
with prompt J/ip production. The muon and kaon tracks are constrained to come from a 
common point of origin and the mass of the fi~^ fi" pair is constrained to the known J/i/j 
mass. Since the intrinsic width of the J/ip is significantly smaller than our experimental 
resolution, the mass constraint improves the resolution of the reconstructed B mass. 

The Pt of each B candidate is required to be greater than 6.0 GeV/c. The proper decay 
length is required to be greater than 100 fim to suppress backgrounds associated with prompt 
J/ip mesons. The signed proper decay length in the B rest frame is defined as 



,.f(m _ ^J/^-Pt 1 _ MbXj/^ ■ Pt , . 



where 



Xj/rp = {xj/^ - xpv)i + {yj/i, - ypv)j (2) 

and [P'-f)^ is the relativistic boost of the B meson. The {xj/^,yj/^) are the transverse 
coordinates of the J/ip decay vertex, and the {xpv,ypv) are the transverse coordinates of 
the event primary vertex. The intersection of the muon tracks as measured in the SVX 
determines the location of the B meson decay. 

The B^ candidate mass distribution is shown in Figure |l|. The distribution is fit with a 
Gaussian signal function plus a linear background using an unbinned maximum likelihood fit. 
The region below 5.15 GeV/c^ has been excluded from the fit since it includes contributions 
from partially reconstructed higher-multiplicity 5-decay modes. The fit yields 387 ± 32 B^ 
mesons. 

V. DIFFERENTIAL CROSS SECTION 

To measure the differential cross section, we divide the B candidate sample into four p^ 
ranges. The invariant mass distributions for each of the pt ranges are then fitted using an 
unbinned maximum likelihood fit which is described in Section A. The determination of the 
geometric acceptance, the efficiencies and the luminosity are described in Sections B, C and 
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FIG. 1: invariant mass distribution reconstructed from the decay — > J/ij^K^. The curve 
is a binned fit to a Gaussian distribution plus linear background and is for illustration only. 

D respectively. The systematic uncertainties are discussed in Section E, and the results are 
presented in Section F. 



A. Fitting Technique 

To measure the meson differential cross section as a function of px-, the B candidate 
sample is divided into four px bins: 6-9, 9-12, 12-15, and 15-25 GeV/c The invariant mass 
distribution for each of the px ranges is then fitted using an unbinned maximum likelihood 
fit to determine the number of B candidates in each px range, as shown in Figure |^. The 
likelihood function is a Gaussian signal plus a hnear background: 

= ~jr^ Jsig H J back 

^^total ^^total 
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where the free parameter Nsig is the number of signal events and Ntotai is the total number of 
candidates in each momentum bin. The function fsig is the Gaussian signal mass function: 



2 



1 ILlzE. 

/,,, = ^^e"n J (4) 



2Tcsai 

where Mj is the candidate mass obtained from a kinematic fit of the muon and kaon tracks. 
The uncertainty CTj on the mass is scaled by a free parameter s in the unbinned maximum 
likelihood fit which is typically ^ 1.2. The parameter M is the mean B mass obtained by 
fitting Figure |I|. The background mass function is linear: 

hack = ^ (^^ - f ) + ^ (5) 

where b is the slope of the background and w is the mass range in the fit (5.15 to 6.0 
GeV/c^). The region well above the B mass yields a better estimate of the slope of the 
background since it is not affected by partially reconstructed B decays. The likelihood 
function is minimized with respect to the parameters Ngig, s and b. The fit yields 160 ± 23, 
114 ± 17, 62 ± 13 and 71 ± 10 events in the four transverse momentum bins. 

B. Acceptances and Trigger Efficiencies 

The acceptance is determined from a Monte Carlo simulation based on a next-to-leading- 



order QCD calculation |l|] using the MRST parton distribution functions [0. The 6-quark 
pole mass is taken to be 4.75 GeV/c^. The b quarks are produced in the rapidity range 



\yb\ < 1.1 with prib) > 5.5 GeV/c. The renormalization scale is /i = /io = yrnf -l-pf'(fc), 
and the fragmentation scale is equal to the renormalization scale. The fragmentation into 



B mesons is modeled using the Peterson fragmentation function with the parameter 



ep set to 0.006 |]T6[. This value was extracted in a fit to data collected at e~^e colliders. 



Recent results from LEP and SLD suggest that lower values of ep and other functions better 



describe the fragmentation of b quarks into B hadrons [|13]- Futhermore the assumption 
that a fragmentation function extracted from e^e^ data is an accurate description of b 
fragmentation at a.pp collider lacks a strong theoretical basis 0. However, the uncertainties 
due to these factors are expected to be smaller than the uncertainty on the renormalization 
scale. 

Decays of Monte- Carlo- generated B mesons into the J/ip and kaon final states are per- 



formed using a modified version of the CLEO Monte Carlo program [|18| which accounts for 

14 



160 - 



— I 1 1 1 1 1 1 p 

Nr = 160 ± 23 



9-12 GeV/c 




150 



— I 1 1 1 1 1 1 1— 

N« = 1 1 4 ± 1 7 



5.80 



B Candidate Mass (GeV/c ) 



B Candidate Mass (GeV/c ) 



o 

> 



-1 1 1 1 1 1 [— 

Nr = 62 ± 13 




J I I I I I I L 



B Candidate Mass (GeV/c ) 



B Candidate Mass (GeV/c 



FIG. 2: candidate mass distribution for the four px ranges. The curve is a binned fit to a 
Gaussian distribution plus hnear background and is for illustration only. 

the expected J/ip longitudinal polarization. Once the B mesons are generated and decayed 
into their final state, a simulation of the CDF detector is utilized. A simulation of the 
trigger efficiency has also been included in the acceptance calculation. The events are then 
processed by the same analysis code used on the data to determine the combined acceptance 
and trigger efficiency for each momentum bin. The Run lA and IB results which incorporate 
different trigger requirements are listed in Table | together with the combined results. The 
uncertainties given are statistical only. 
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TABLE I: The product of the trigger efficiency and the acceptance in the px bins for Run lA, Run 
IB and the integrated luminosity- weighted average for Run 1. 



Pt range 


Trigger efficiency x acceptance (%) 


(GeV/c) 


Run lA 


Run IB 


Run 1 


6-9 


2.01 ± 0.02 


1.61 ± 0.02 


1.70 ± 0.02 


9-12 


5.29 ± 0.05 


4.20 ± 0.04 


4.44 ± 0.03 


12-15 


8.36 ± 0.10 


6.53 ± 0.09 


6.93 ± 0.07 


15-25 


11.96 ± 0.14 


9.26 ± 0.12 


9.86 ± 0.10 



C. Efficiencies of the Additional Selection Requirements 

The detector acceptance and trigger efficiencies described in the previous section did not 
account for all of the criteria for selecting a B candidate. The efficiencies of the additional 
selection requirements are discussed in this section. Most of these efficiencies are determined 
using large CDF data samples. 

There are two components that comprise the tracking efficiencies. The first part is the 
efficiency of the tracking in the Level 3 trigger system which is determined using an inclusive 
single muon data set. The efficiency is measured to be (97 ± 2)% for Run IB. During Run 
lA, a portion of the data-taking suffered from the start time of each event being incorrectly 
determined. The result was an inefficiency in reconstruction at Level 3 which was determined 
to be ~ 4% [|l^ averaged over all of Run lA. The Level 3 Run lA efficiency is (93 ± 2)%. 

Once an event has been accepted at Level 3, one must account for the offiine CTC track 
reconstruction which may improve the muon track quality or find new tracks that are missed 
at Level 3. It is also necessary to correct for the track finding efficiency for the kaon track 
since it is not required in the Level 3 trigger. A detailed study pO[ of the CTC track 



reconstruction efficiencies was conducted. To measure the efficiency, we simulate single kaon 
tracks with the CDF Monte Carlo. We then combine the generated CTC hits for such a 
kaon with the hits in an event with an identified displaced J/ ip from the CDF data sample. 
Hits in the CTC are characterized by a leading edge and a time-over-threshold. Where a real 
and simulated hit overlap, the hits are combined. Thus the leading edges used in the track 



16 



reconstruction may be obscured for the simulated kaon as they would be for real particles. 
We then run the full track reconstruction program on the modified event and search for a 
track corresponding to the embedded kaon. We find the efficiency of the track reconstruction 
to be (99.61^09)% for particles with pt > 0.8 GeV/c that traversed all layers of the CTC, 
independent of instantaneous luminosity. The Run lA single track reconstruction efficiency 
of (98.5 ± 1.4)% is taken from Ref. §. 

The muon segment reconstruction efficiency is found to be (98.0 ± 1.0)% resulting in a 
combined efficiency of (96.0 ±1.4)%. The efficiency of requiring both muons from the J/ip to 
have a muon chamber track segment that matches a track reconstructed in the CTC is found 
to be (98.7 ±0.2)%. The efficiency of this cut is determined from a sample of J/ip candidate 
events containing muons that were required to pass less stringent matching requirements at 
Level 3. 

The fraction of events in which both muons from the J/ip have been reconstructed in the 
SVX is measured using a large J/ip data set. This fraction is (52.4 ± 0.6)% for Run lA and 
(56.3 ± 0.2)% for Run IB. The fraction for Run IB is larger than Run lA because the inner 
layer of the SVX detector was moved closer to the beamline, eliminating a small separation 
between silicon wafers in the first layer present in Run lA. 

The efficiency to reconstruct a B meson with a proper decay length ct greater than 100 
/xm is determined using Monte Carlo simulations. The ct resolution is measured in the 
J/ip data set by fitting the proper lifetime of events in the sidebands of the B candidate 
mass distribution with a Gaussian function for the prompt component and an exponential 
function for the long-lived component. The lifetimes of the Monte Carlo generated events 
are then smeared using the resolution measured in each px range. The efficiency showed no 
significant variation with the B transverse momentum even though the proper ct resolution 
was degraded by a factor of 2 from the lowest to the highest pt bin. The efficiency of 
(78.4 ± 0.5)% is the mean of the values measured in each pt bin. 

The reconstruction efficiencies are summarized in Table For the B candidates de- 
caying to particles completely contained within the detector acceptance, the reconstruction 
efficiency is (36.4 ± 1.0)%. 
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Source 


Efficiency in % 




RunlA 


Run IB 


CTC tracking 


(98.5 ± 1.4)3 


(99.6l°i)3 




= 95.6 ± 2.1 


= 98.8;?:^ 


L3 /x"*"//" tracking 


93 ±2 


97 ±2 


CTC-M linking 


(99.8 ± 0.2)2 
= 99.6 ±0.3 


Muon chamber 


(98.0 ± 1.0)2 


efficiency 


= 96.0 ± 1.4 


fj,^fj-^ matching cut 


98.7 ±0.2 


Z vertex cut 


95.3 ± 1.1 


93.7 ± 1.1 


SVX fraction 


52.4 ±0.6 


56.3 ±0.2 


(i > 100 //m 


78.1 ±0.5 


Total 


36.4 ± 1.2 



TABLE II: Summary of reconstruction efficiencies for the B meson. The efficiencies that are not 
common between lA and IB are averaged and weighted by integrated luminosity. 

D. Luminosity Determination 

At CDF the luminosity is measured using two telescopes of beam-beam counters to an 
accuracy of about 4%. We studied the quality of the integrated luminosity calculation in the 
inclusive J/ip ^ sample. After correcting for the time-dependent trigger efficiency, 

we found that in Run IB the measured J/ip cross section fell linearly as a function of 
instantaneous luminosity C However, for any narrow range of JC, was constant as a 
function of time. Since the minimum luminosity of the data sample is 4 x lO^*^ cm^^s"^, we 
have considered two possible extrapolations of as a function of £ to £ = to calculate 
a corrected integrated luminosity. The first extrapolation is performed assuming that the 
linear dependence is valid below £ < 4 x 10^° and that: 

JC'dt^ Jm^^dt (6) 
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We also perform the extrapolation assuming that no correction is needed below £ < 4 x 10^°. 
The luminosity correction is taken to be the average of the two extrapolations and we assign 
a systematic uncertainty that covers the range between the two hypotheses. The correction 
to the integrated luminosity for Run IB is 

Rc = I Cdt/ J C'dt^ 0.88 ±0M. (7) 
E. Systematic Uncertainties 

We divide the systematic uncertainties in the measurement of the meson production 
cross section into two classes: px dependent uncertainties {systpj.) that change from one pt 
bin to the next and fully correlated uncertainties (systfc) that are independent of pr- 

1. Pt dependent systematic uncertainties 

The px dependent systematic uncertainties include variations of the production and decay 
kinematics that would affect the determination of the acceptance. We have considered effects 
due to the model used to generate the b quark spectrum and uncertainties in our knowledge 
of the trigger efficiency. 

The model used to generate the b quarks is based on a QCD calculation at next-to- 
leading order. Large uncertainties in the calculation are due to unknown higher-order effects. 
These effects are quantified by estimating the scale dependence when the renormalization 
and factorization scales are varied by a factor of 2 above and below their central value of 
II — fj,o — ^Jp^ -\- ml. The Peterson fragmentation parameter is varied by ±0.002 around 
its central value of ep — 0.006. In each case the uncertainty on the acceptance is taken to 
be the difference between the acceptance found with the central value and the value found 
when each variable is varied by the indicated amounts. The dependence of acceptance on 
the parton density parametrization and the b quark mass are much smaller and are not 
included in the systematic uncertainty. In addition, the parameters of the trigger simulation 
are varied by ilcr. The total pt dependent uncertainty is given by the sum in quadrature 
of the Pt dependent systematic uncertainties summarized in Table III. 
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2. Correlated systematic uncertainties 



The correlated systematic uncertainties include uncertainties that are independent of the 
B meson px spectrum. The largest of these uncertainties is due to limited knowledge of 
the — s> J/ipK^ branching ratio which yields a systematic uncertainty of about 10%. 
Other sources of correlated uncertainties are due to the uncertainty on the total reconstruc- 
tion efficiency shown in Table II and knowledge of the integrated luminosity collected at 
CDF during Run 1. There is an additional systematic uncertainty associated with the re- 
construction of kaons that decay inside the CTC volume. A simulation shows that about 
8% of the kaons decay in flight, of which half are successfully reconstructed [Q. We assign 
the full value of the correction as an uncertainty for the kaon acceptance of (96 ±4)%. This 
assumes that such tracks are modeled realistically in the simulation. The total correlated 
uncertainty of 12.7% is given by the sum in quadrature of the fully correlated systematic 
uncertainties summarized in Table IV. 



TABLE III: Summary of pT dependent systematic uncertainties. 



Source 


Fractional uncertainty in each px bin 


Pt range (GeV/c) 


6-9 


9-12 


12-15 


15-25 


QCD renormalization uncertainty 


1.6% 


1.5% 


1.7% 


1.5% 


Peterson parameter uncertainty 


0.7% 


1.6% 


1.0% 


1.7% 


Trigger efficiency uncertainty 


3.1% 


2.7% 


2.1% 


1.7% 


Pt dependent total (syst^j,) 


3.6% 


3.5% 


2.9% 


2.8% 



F. Results 

The differential cross section da/dpr is calculated using the following equation: 

dpT ApT-C'-A-e-B ^ ^ 

where Ngig is the number of charged B mesons determined from the likelihood fit of the mass 
distribution in each p^ range. The factor of 1/2 is included because both B~^ and B~ mesons 
are detected while we report the cross section for B^ mesons assuming charge invariance in 
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TABLE IV: Summary of fully correlated systematic uncertainties. 



Source 


Fractional uncertainty 


Reconstruction efficiency 


±2.7% 


Luminosity uncertainty 


±4.1% 


Luminosity correction 


±4.5% 


Branching ratio uncertainty 


±10.2% 


Kaon decay-in-flight uncertainty 


±4.0% 


Fully correlated total (syst/c) 


±12.8% 



the production process. The width of the pt bin is Ap^ and C is the corrected integrated 
luminosity of the sample. The geometric and kinematic acceptance A is determined from the 
Monte Carlo simulation and includes the kinematic and trigger efficiencies. The efficiency 
e is the additional reconstruction efficiency not included in the simulation. The product of 



branching ratios B is determined using the the world-average branching fractions: 



BR{B^ J/^ K^) = (10.0 ± 1.0) X 10"^ (9) 

BR{J/ip /iV") = (5.88 ± 0.10) X 10"^ (10) 

Table V lists the differential cross section as a function of px- The three uncertainties 
quoted on the cross section are statistical (stat), px dependent systematic (systp^), and 
fully correlated systematic (syst/c), respectively. 



TAB] 


V: meson differential cross section from the Run 1 data. 


(pt) 
(GeV/c ) 


Events 


Acceptance 

(%) 


Cross section 
(nb/[GeV/c]) 


7.34 


160 ± 23 


1.70 ± 0.02 


815 ± 117(stat) ± 31(systp7,) ± 104(syst/c) 


10.35 


114 ± 17 


4.44 ± 0.03 


222 ± 33(stat) ± 8(systp7,) ± 28(syst/c) 


13.36 


62 ± 13 


6.93 ± 0.07 


77.5 ± 16.2(stat) ± 2.4(systp^) ± 9.9(syst/c) 


18.87 


71 ± 10 


9.86 ± 0.10 


18.7 ± 2.6(stat) ± 0.6(systpy) ± 2.4(syst/c) 
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Figure ^ shows the measured differential cross section at the mean px of each bin com- 
pared to the next-to-leading-order QCD [|1| calculation using the MRST parton density 



functions |]T4|. The experimental points are plotted at (pr) which is the value of pr for 
which the theoretical differential cross section equals the mean cross section in each 
momentum range 



da 



^PT da , , , 

dpT (11) 



dp' 

The dashed lines in Figure ^ indicate the change in the theoretical predictions as the 
b quark mass is varied between 4.5 and 5.0 GeV/c^, the renormalization scale is varied 
between /io/2 and 2/io, and the Peterson fragmentation parameter is varied between 0.004 
and 0.008. The solid curve is for the central values of these parameters: rrib = 4.75 GeV/c^, 



Ho = v^^+Pt, and ep = 0.006. The fraction of b quarks that fragment into is 



fu = 0.375 ± 0.023 pg. This fraction is varied between 0.352 and 0.398. 

The comparison between data and theory for da/dpx is aided by plotting the ratio of 
data/theory on a linear scale, as shown in Figure ^. The level of agreement between the data 
and the theoretical prediction is determined by fitting a line through the four ratio points. 
The fit yields a scale factor for data/theory of 2.9 ± 0.2 (stat © systp^) ± 0.4 (systfc) with a 
confidence level of 72%. The first uncertainty on the scale factor is the uncertainty returned 
by the fit to the ratio points whose uncertainties were determined by summing the statistical 
and the px dependent systematic uncertainties in quadrature. The second uncertainty is the 
fully correlated systematic uncertainty. The hatched band shows the magnitude of the fully 
correlated uncertainty which arises mainly due to the poor knowledge of the J / ipK^ 

branching franction. Also shown is a comparison between the shape of the QCD predictions 
obtained using a different set of parton distribution functions determined by the CTEQ 
collaboration |2^. The effect of changing the parton distribution functions is negligible 
in comparison with the variation associated with uncertainties in the b quark mass, the 
fragmentation parameter and the renormalization scale shown by the dashed curves. 



VI. THE TOTAL CROSS SECTION 



The total cross section is obtained by using a method similar to the one used for the 
determination of the differential cross section. However, the last tranverse momentum bin, 
15-25 GeV/c, is replaced with the invariant mass distribution for B"^ candidate events with 
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FIG. 3: -B"*" meson differential cross measurements compared to the theoretical prediction. The 
solid curve is the theoretical prediction for nih = 4.75 GeV/c^, /xq = y^m^ + Pt-: = 0.006 and 
fu = 0.375. The dashed lines illustrate the changes in the theory once these parameters are varied 
as explained in the text. 

Pt > 15 GeV/c shown in Fig. |[ With 81 ± 11 candidates and an acceptance of (10.19 ± 
0.16)%, the integrated cross section for px > 15 GeV/c is 207 ± 28(stat) ± 5(systpy) ± 
26(systjc) nb. The integrated cross section for B transverse momentum px > 6.0 GeV/c 
and \y\ < 1.0 is given by: 

a(B+) = y ^^/^ (12) 

where Ni is the number of charged B candidate events in each momentum bin, Ai is the 
acceptance and e is reconstruction efficiency. The total cross section is: 

(Tb{pt > 6.0 GeV/c, \y\ < 1.0) = 3.6 ± 0.4 (stat © systp^) ± 0.4 (systfc) /xb. (13) 

where the first uncertainty is the sum in quadrature of the statistical and pt dependent 
systematic uncertainty, and the second uncertainty is the fully correlated systematic uncer- 
tainty. 
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FIG. 4: Plot of data/theory as a shape comparison with the NLO QCD differential cross section 
calculations. 

VII. SUMMARY 



The exclusive decay ^ J/ip has been used to measure the production cross section 
of the meson from data collected by the CDF detector. A sample size of 387 ± 32 events 
is obtained from / Cdt = 98 ± 4 pb^^ of 1.8 TeV pp collisions produced by the Fermilab 
Tevatron collider. 

The measured total B^ production cross section for Pt{B) > 6.0 GeV/c and \y\ < 1.0 is 



asiPT > 6.0 GeV/c, \y\ < 1.0) = 3.6 ± 0.6(stat © syst) fih 



(14) 



where the uncertainty is the sum in quadrature of the statistical and both correlated and 
Pt dependent systematic uncertainties. The differential cross section is measured to be 
2.9 ± 0.2 (stat © sjstprp) ± 0.4 (systfc) times higher than the NLO QCD predictions with 
agreement in shape. The first uncertainty is the sum in quadrature of the statistical and pt 
dependent systematic uncertainty and the second is the correlated systematic uncertainty. 
The new measurement of the B~^ differential cross section confirms that the absolute rate is 
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FIG. 5: candidate mass distribution for pxiB) > 15 GeV/c. The curve is a binned fit to a 
Gaussian distribution plus linear background and is for illustration only. 



larger than the hmits of that predicted by typical variations in the theoretical parameters. 
These measurements supersede those of reference 0. 
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